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Abstract 

Given a gapped Hamiltonian of a spin chain, we give a polynomial-time algorithm 
for finding the degenerate ground space projector. The output is an orthonormal 
set of matrix product states that approximate the true ground space projector up to 
an inverse polynomial error in any Schatten norm, with a runtime exponential in the 
degeneracy. Our algorithm is an extension of the recent algorithm of Landau, Vazirani, 
and Vidick [1] for the nondegenerate case, and it includes the recent improvements 
due to Huang [2]. The main new idea is to incorporate the local distinguishability of 
ground states on the half-chain to ensure that the algorithm returns a complete set of 
global ground states. 

Strongly correlated quantum systems are at the heart of such diverse physical phenom¬ 
ena as frustrated magnets, high-temperature superconductors, and topological quantum 
phases, to name but a few, and understanding their low-temperature behaviour is a grand 
challenge at the interface of physics and computer science. The computational tools to 
study these systems numerically are however inevitably stymied by the exponential growth 
of the Hilbert space of an n-particle system. 

Many numerical tools have been developed to surmount the computational difficul¬ 
ties surrounding the exponential growth of the state space. Tools like quantum Monte 
Carlo [3,4], coupled cluster methods [5], and the density-matrix renormalisation group 
(DMRG) [6] all exploit some underlying physical structure to perform heuristically effi¬ 
cient simulations of interacting quantum systems. However, like all heuristics they often 
lack basic theoretical guarantees except in certain special cases, and some of them have 
failure modes that are difficult to diagnose or characterise. This makes it hard to precisely 
determine which physical features enable efficient simulation and which features might be 
responsible in the cases where the heuristic fails. 

Very recently, a success story has emerged in the rigorous understanding of gapped 
ID spin chains. Hastings’ proof of the ID area law [7] established that approximating the 
ground state of a gapped spin chain is in NP. The proof technique crucially used the idea of 
matrix product states (MPS) [8], an ansatz for quantum systems that exploits the area law 
structure of entanglement [9]. The hrst rigorous nontrivial algorithm for hnding ground 
states, due independently to Schuch and Cirac [10] and Aharonov, Arad, and Irani [11], 
used dynamic programming to hnd MPS ground states. The runtime of the algorithm was 
exponential in the bond dimension, and since Hastings’ result used a polynomial bond 
dimension to achieve a sufficiently accurate approximation, the runtime could only be 
guaranteed to be exponential for general gapped ID systems. The hrst subexponential¬ 
time algorithm was found by Arad, Kitaev, Landau, and Vazirani [12]. They used the 
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same algorithm as Refs. [10,11], but showed that the ground state could be approximated 
to sufficient accuracy with an MRS of bond dimension where n is the length of the 

chain. This line of research culminated in the breakthrough algorithm of Landau, Vazirani, 
and Vidick (LVV) [1], who gave the first polynomial-time algorithm for approximating 
the ground state in the case of a non-degenerate ground space. This was subsequently 
improved by Huang [2] who gave a deterministic algorithm with an exponentially improved 
scaling with respect to the gap that also applied to highly frustrated systems. 

In this paper, we extend the LVV algorithm and its subsequent improvements by 
Huang to the case where the Hamiltonian has a degenerate ground space. To give a 
precise statement of our results, we first introduce some notation. 

Consider a ID chain of d-level quantum systems (qudits) with d = 0(1) and a Hamil¬ 
tonian H that is a sum of local terms on contiguous sites, and let n be the number of 
sites in the chain. After sufficient coarse-graining, shifting the energy, and dividing by the 
largest operator norm from the terms in the Hamiltonian, we can take the Hamiltonian 
to be of the form H = where 0 < Hi < 1 and Hi acts non-trivially only on 

the i and i + 1 qudits. We additionally require that this transformed Hamiltonian has a 
constant lower bound ei — cq > e > D(l) on the spectral gap, i.e. the difference between 
the first excited energy (ei) and the ground state energy (cq) of H is bounded from below 
by e independent of n. Our Hamiltonians will have a degenerate spectrum in the ground 
space with a degeneracy g = 0(1). We say that such a Hamiltonian is in standard form, 
though this is not a unique normal form because of the ambiguity in the coarse-graining 
step. 

Our main result is an approximation of the ground space in the form of an orthonormal 
basis of matrix product states that approximate the true ground space projector. More 
precisely, we have the following. 

Theorem 1 (Main result). Given an n-qudit Hamiltonian H in standard form and con¬ 
stant bounds e and g on the gap and degeneracy respectively, then for every sufficiently 
small approximation error g satisfying Ifg = Algorithm 1 constructs the orthonor¬ 

mal approximate ground states {\'yj)}j in matrix product form for 1 < j < g in runtime 
, The approximation error gives a bound in the Frobenius norm of 

\\G-T\\^<g, 

where G is the true ground space projector, and T := l7i)(7il approximate 

ground space projector. 

The proof of this statement is given in Section 2.7, and a more precise dependence of 
the runtime on the gap e and the ground space degeneracy g is made explicit in Section 2.8. 
Note that since both G and T have rank g, we can get an approximation error in the trace 
norm at the cost of a factor of ^/^ = 0(1) increase in error. 

Our result relies heavily on the ideas in LVV [1] and the improvements of Huang [2], 
so let us recall briefly how this algorithm proceeds. The LVV algorithm proceeds by 
constructing a viable set, which roughly speaking is a small set of matrix product states 
on a half-chain that is guaranteed to contain the local Schmidt vectors of an approximate 
ground state on the full chain. In each step of the algorithm, the half-chain of the current 
viable set is extended by one site to create a new viable set. This initial extension is too 
large, but trimming away some states with too high an energy yields another suitably 
small set that is still viable. 

The main new idea is an extension of the size trimming step to degenerate systems, 
based on the concept of local distinguishability of ground states. The algorithm works by 
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iteratively sweeping across the chain, and finds the degenerate ground states one-by-one 
in a manner similar to that of LVV on each iteration. As we run these iterations, we must 
ensure that the viable set of potential ground states on the half-chain contains a witness 
that not only has low energy, but also does not simply reduce to something in the span of 
the previously found ground states. To ensure this, we break the analysis into two cases 
of high- and /ow-distinguishability. 

In the low-distinguishability case, the Schmidt vectors on the left half-chain of the 
previously found ground states and the desired witness have large overlap. By recycling 
these Schmidt vectors, we therefore get a viable set for this new witness for free. 

In the high-distinguishability case however, this recycling would not provide a viable 
set. In this case we need to trim the viable set in a way that not only ensures the existence 
of a low-energy witness, but a witness which has low overlap with the existing states. If the 
relevant Schmidt vectors of the desired witness and the existing states have low overlap, 
then we can optimise for a state of low energy and overlap efficiently. 

1 Preliminaries 

Let us begin by giving a more precise statement of the LVV algorithm, with the definitions 
suitably extended to our treatment of the degenerate case. 

Given the Hamiltonian H = X]j=i often consider a partition of the system 

into the first i spins and the rest, we will call this the z|z -|- I cut. The Hamiltonian on 
a contiguous subset of spins {j,..., /c} is denoted = Yli=j notation the 

local terms can also be expressed as Hi = and the full Hamiltonian as H = 

In the special case that an index i is clear from context, we define Hi = and 

Hr = , 1 ]. Note that this implies that H = Hr + Hi + Hr. The Hilbert space on 

which acts is denoted by Hy^k]- Similarly, we have the Hilbert space H = Hy^n]-, 

similarly for the other possible subscripts. It will sometimes be convenient to subtract a 
constant energy shift from a Hamiltonian, in which case we define H'j^ = Hr — ei, where 
cl is the least eigenvalue of Hr, and similarly for H'j^ and eR. When a particular state 
p is clear from context, we refer to Tii[pHR) as the left energy of p, and similarly for the 
right energy. 

A viable set is simply a set of vectors that can be used to construct the Schmidt vectors 
of a new approximate ground state on a half-chain, assuming we already have approximate 
ground states I 71 ),... ,\'yh-i)- The aim is to inductively grow the region of viability to the 
entire system starting from one end of the chain. The span of such a set gives the desired 
small domain on which candidate ground states can be found efficiently. 

Definition 2 . (Viable Set) A set of states S is {h,i,s,b,6 or A)-viable if: 

• The states in S are supported on the first i < n qudits of the system, i.e. S C 

• The cardinality is bounded |5| < s. 

• Each state in S is represented by an MPS with bond dimension at most b. 

• There exists a witness state |u) which has, depending on context, either 

— {5 case) overlap with the ground space projector G satisfying ||G|u)|| >1 — 5, 
or 

— (A case) the expected energy satisfying {v\H\v) < cq -|- Ae. 

• The reduced density operator of this witness must be supported on Span(5). 
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• The witnesses state lu) is orthogonal to all previously constructed approximate 
ground states I 71 ),... ,\'yh-i)- 

When it is clear from the context (and it almost always will be clear), we will drop 
the explicit h dependence in the definition of viable set and refer to a {i, s, b, <5)-viable set, 
and similarly for A. 

The two error parameters specified, 5 or A, will be used contextually and we will 
sometimes refer to them as the overlap error and the energy error, respectively. The 
two error parameters however can be used interchangeably, though at the cost of some 
overhead which we would like to avoid. The proof is simple and can be found in [1,2]. 

Lemma 3 (Interchangeability of 6 and A [1,2]). For the ground state projector G of a 
Hamiltonian H, and any |u) ^H, we have 

{v\H\v) < eo + Ae ||G |u) || > 1 — A 

||G |u) II > 1 — h {v\H\v) < eo + 25n 


which allows us to convert between the two measures of viability 


A)-viable (•, •, 6 = A)-viable 

2n 

6)-viable (•,•,•,•, A = — 5)-viable. 

For later use, we define the following absolute constants c and d that are independent 
of all other parameters, and a simple function / of the ground state degeneracy, 

c :=6.25 X 10"^^ c' :=10"i^ / :=g{2g + 1), 


where we note that 200-^ = c'. These values can be used to make all of the constants 
hidden by the big-0 notation in our derivations and procedures completely explicit. How¬ 
ever, our final result still relies on certain implicit estimates (in particular, the estimates 
in Refs. [13,14]) so we are not able to give a completely explicit scaling. 


2 The Algorithm 

Algorithm 1 works by first running the slight modifications on the existing nondegenerate 
ground state approximation algorithms of Ref. [1,2], which is discussed in more detail in 
Section 3.1.1. As discussed earlier, the extended algorithm works by iteratively finding 
ground states one-by-one, with each guaranteed to be orthogonal to those previously found 
until all g ground states have been found. It is not necessary to know a priori a specific 
value for g, but unless g = 0(1) holds, then there is no guarantee from the area law 
alone that all the ground states can be approximated well by the algorithm. Moreover, if 
the algorithm is iterated more times than the true degeneracy of the system, then it will 
necessarily produce high energy states which will herald this over-estimation, as we will 
show in Appendix A. As such we will henceforth assume that the degeneracy g is known. 

The Step/FinalStep procedures will be discussed below in Section 2.1, and its four 
sub-procedures will follow in Sections 2.2-2.5. Finally the ApproxGroundState proce¬ 
dure will be discussed in Section 2.6. 
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Algorithm 1 Ground space approximation 
1: procedure DEGENERATEGSA(ff, g, e, rj) 

2: | 7 i) ^ NONDEGENERATEGSA(ff, e, rj) t> Run non-degenerate algorithm 

3: i? •(— ( 7 i|i/| 7 i) 0 Approximate the ground energy 

4: for all h £ [2,5'] do 

5: Shft {1} Initial viable set 

6: for all i G [1, n — 1] do 

7: Ei £- ^ 7 i j] | 7 i) > Estimate energy of half-chain 

8: SGHMiDTVECS(i, | 7 i),.. .,|7h-i)) t> Extract left Schmidt vectors 

9: Sh,i ^ STEP(5/i^j_i, Lh^i, H, e, E, Ei) t> Step viable set one qudit along 

10: end for 

11: Sh^n ^ FinalStep(5/i^„_i, H, e, E, rj) > Low-error final step 

12: | 7 /j) •(— AppROxGROUNDSTATE(Lf, Sh^ni Iti): • • • j l7h-i)) » Extract witness 

13: end for 

14: return jyi) , ... , 

15: end procedure 


Algorithm 2 Stepping functions 
1: procedure STEP{Sh,i-i, Lh,i, H, e, E, Ei) 
2: ^ EXTEND(S’/j,i_i) 

3: ^ Trim(4;J, Lh,i, H, e, Ei) 

^ Trungate(S’[^], e) 

5: £- REDUGE(5g, H, €, E) 

6: return 5^"^] 

7: end procedure 


0 Increment i 
> Reduce cardinality of set 
i> Truncate state bond dimensions 
0 Push back down the error level 


i: procedure FinalStep(S'/j^„_i, ff, e, E, g) 


9: 

o(l) , 

- EXTEND(S’h,n-l) 


> Extend to final qudit 

10: 

o(2) , 
‘^h,n 

- Tmuisi^l, Lh,n, H, 

C) En) 


11: 

0 ( 3 ) , 
‘^h,n ^ 

- Trungate(S'[^^, e) 



12: 

o(4) , 
‘^h,n ^ 

- FinalReduge(S'[^^, 

H, e, E, g) 

0 Reduce error to desired level 

13: 

return 




14: end procedure 


2.1 Step 

The Step procedure, shown above in Algorithm 2, takes a viable set defined on i — 1 qudits, 
and constructs a new set which is viable on i qudits. As well as physically extending the set, 
this procedure can be efficiently performed such that neither the cardinality nor the error 
grow. Starting from the trivially (0,l,l,0)-viable set {!}, the Step procedure allows you 
to inductively construct (i,s,6,A = ce®//'^)-viable sets Sh,i, for all i. Here the parameters 
s and b are all polynomial functions of n and are given in Table 1. On the final qudit the 
Step procedure is tweaked in the form of FinalStep, which performs a more powerful 
error reduction procedure to bring the error down to the desired final level. 

By constructing a viable set on the entire system, a polynomial-sized subspace has 
been constructed on which the Hamiltonian may be efficiently optimised, allowing the 
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Table 1: The viability parameters of the sets produced during Lines 8-14 of Algorithm 1. 
Red indicates the desired parameter change in each sub-step, grey the weaker of the two 
error parameters as given by Lemma 3. Here p is some absolute polynomial p = and 
the remaining polynomials scale as po,pi,p 2 ,q = where a = 0(l + e“^y^log77“^/log n). 
The polynomials s := ppi + pq + q and b := pp 2 denote the fixed points of the cardinality 
and bond dimension under the Step procedure. 


previously inaccessible witness state to be harvested, yielding the hnal desired ground 
state approximation. This hnal optimisation is performed in ApproxGroundState, 
which follows the SxEP-based viable set construction. 

2.2 Extension 

By induction, suppose we start with a (i — l,s,b,A = ce®//^)-viable set Si-i. The hrst 
step is to take a viable set on z — 1 qudits and extend it trivially to a viable set on i 
qudits. To do this we take the element-wise tensor product of the existing viable set with 
a basis on the zth qudit. It is trivial to see that the set resulting from Extension, 
is {i,ds,b,A = ce®//^)-viable; see Ref. [1] for the explicit argument. We simply note that 
the net result of this step is to increase the cardinality of the viable set from s to ds, while 
keeping all other parameters hxed. 

2.3 Size trimming 

For the size-trimming we are going to use two nets to compensate for our ignorance of 
two properties of the ground states: the local energy differences between them, and their 
entanglement structure across the i\i -I- 1 cut. The former is represented by the difference 
in the expectation value ot between our desired witness and (say) the first ground state 
found. The latter is represented by the boundary contraction (dehned below) of the witness 
state in question. In both cases we are going to construct a net on the space of possible 
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values of these unknowns, and perform the necessary optimisations in the neighbourhoods 
of these elements, thus allowing for error guarantees in spite of this ignorance. 

2.3.1 Boundary contraction 

Viable sets, like reduced density matrices, capture local information about a given state. 
As defined, viable sets are solely determined by the left Schmidt vectors of a witness 
state; the reduced density operator however also contains the Schmidt coefficients across 
the cut. As such, merely finding vectors of low energy on the half-chain [l,i] in no way 
guarantees consistency with a low-energy global state, as they might not be consistent 
with the entanglement structure across the i\i + 1 cut. As well as this the system may 
be frustrated, with a globally low-energy state appearing locally excited. Both of these 
problems are circumvented by performing the desired optimisations alongside boundary 
contraction constraints. 

Given a state of Schmidt rank B and Schmidt decomposition across the i\i + 1 cut 
given by |u) = l®i) \^j)i —)■ T-Lr be the partial isometry specihed by 

Uv \j) = \bj)- Then we have the following. 

Definition 4 (Boundary Contraction [1]). The left state of |u) with respect to the i\i + 1 
cut is defined as 

|ls(u)) := Ul |u) = ^ Xj \aj) \j) G (g) . 

3 

The boundary contraction of |u) is dehned as 

cont(u) := Tr[^j_i] (|ls(u)) (ls('(;)|) = C/J Trp_j_i] {\v){v\)U^. 

Note that the density matrix cont(u) is supported on T-Li (g) C^. 

The interpretation of a boundary contraction is the following. The difference in right 
energies between a candidate state and a valid witness state for the full chain is bounded 
by how close the boundary contractions of the states are. By extending a state by the 
as-yet-unknown witness state, we can argue the existence of a low energy solution within 
our extended viable set. 

Lemma 5 ( [1,2]). Let a be a density matrix on and |u) = Ylf=i l®i) \bj) £ T~L 

be a Schmidt-deeomposed state. The density matrix a' := U^aUl on Ti satisfies 

Tr {a'H) < TT{aHL) + {v\{Hi + Hr)\v) 

+ ||Trp_j_i] (a) - cont(u)||^ (1 -F - {v\Hr\v)) . 

where n„ = UyUv is the rank B projector from the isometry Uy. 

Proof. We give a sketch of the proof and refer the reader to Ref. [1] for more detail. 
Begin with TT{a'H) = Tr[((T' -|- |t')(u| — |u)(u|)^f] for some ground state |u). Expanding 
H = Hr -\- Hi + Hr and using the matrix Holder inequality together with the bound 
\\Hi\\ < 1 yields the proof. □ 

2.3.2 Nets 

Recall that an e-net on a metric space is a set of points such that any point in the space is 
within distance e of a point in the net. We will use two nets in our argument: one on the 
interval [—1,1 -|- ??] in the standard absolute value metric and one on the set of matrices 
with bounded max-norm on C'^ (g) C'® in the trace norm. 
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We denote by Airj the r//2-net on [—1,1 + r/] whose points are 

-! + ?/,••• , 1 - 1,1 + ??} • 

The cardinality of Ain is at most [2/r/] + 1, and for any real value x G [— 1 — ry, 1 + r/] 
there exists a y G Ain such that the following one-sided error bound holds, 0 < y — x < r]. 
We refer to this net as the energy net. The choice of limits [—1,1 + ry] is justified in 
Appendix C. 

Let Afn be the matrices on ( 8 ) C'® with entries whose real and imaginary parts are 
in the set 

{-l,-l + 7y/i?V,--- ,l-ri/B^d\l} 

The cardinality of J\fn is upper bounded by (2 [Bd/ry] + 1)^ . For any positive semidef- 

inite matrix p on C'^ (8 with trace at most 1, there exists a point X G Afn such that 
Up — < Tj. We refer to this net as the boundary contraction net. 

2.3.3 Procedure 

Let Lh,i be the set of left Schmidt vectors with respect to the i\i + 1 cut of the states 
Iti) ) • • •) l 7 h)- Furthermore define to be the projector onto Span(L/i j). Now define 
parameters 5 = 8 y^e^//^, ry = AdeJf, and ^ to be chosen later in Eq. (7). 

Trim: For each boundary contraction X G A /'^/2 with bond dimension (to 
be defined below in Eq. (3)) and left-energy difference Y G Ain^ take cr to be a 

density matrix supported on SpanU ® C^. We will specifically choose 
(T to be a solution to the degenerate size-trimming convex program: 

min Tr(n/i^jCj) ( 1 ) 

where Tr [Hl (cj - I 71 ) ( 71 1)] <Y, 

Tr(cj) = 1, cr > 0. 

Now consider the set of vectors {\vk)}k consisting of all the eigenvectors of a with 
eigenvalues at least 10~^/g. Now define 

~ ’ ( 2 ) 
where \vkj) is the jth left Schmidt vector of \vk) across the i\B cut. 


The next Claim, whose proof is postponed until the next section, establishes that the 
viability parameters after each Trim step remain bounded as advertised. 

Claim 6 (Proof in Section 3.2). The set resulting from Trim, S^^\, is {i,pi,dsb -\- ,5 = 

-viable, where 

Pii'n), q{n) = where a = 0(l -|- ylogry^^VloS^) ■ 

2.4 Bond truncation 

Eirst, we recall some results on low-rank approximate ground states [12,14,15]. Let us 
define a function 

Bs := iog3/4 5-1) _ 
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Specifically we will take Bs to be defined by the following lemma, a corollary to the area 
law [7,12,14], 

Lemma 7 ([12,14]). For any ground state |r) and bipartite cut of the system in question, 
there exists a state [y) with Schmidt rank across that cut bounded by Bs such that |(r| 7 )| > 
1 — 5. Moreover, there exists an MPS approximation ly^) with bond dimension across every 
cut of at most Bs/n and for which |(r|y')| > 1 — 5. 

If we know the Schmidt decomposition of a given state across a bipartite cut, then the 
Eckart-Young theorem [16] guarantees that the truncation of that state to Schmidt-rank 
D is the best approximation with said Schmidt rank. 

Definition 8 (Truncation). Consider a state |u) with Schmidt decomposition |u) = 
l®j) l^j) across the i\i + 1 cut, where any degeneracy in the Schmidt coefficients 
is resolved arbitrarily, and with non-increasing Aj. Given an integer D, define the trunca¬ 
tion of |u) by 

D 

Trunc)^ |u) := Xj \aj) \bj) . 
i=i 

This approximation is optimal in the following sense [16]. For any state |u), and any 
integer D, the truncated vector \v') := TruncD |u) / ||Trunc£) |u)|| satishes (u|u') > |(u|r(;)| 
for any |t(;) of Schmidt rank at most D. 

Combining these results, we get a precise statement that low-energy states can also be 
approximated faithfully with low-rank approximations. See Ref. [2] for a proof. 

Lemma 9 ([2]). For a state |u) with {v\F[\v) < eo-|-Ae where A(l-|-e) < 1/4, the truncated 
state lu)) := Trunc^^^ |u) / ||Trunc^^^ |u)|| has energy {w\H\w) < cq -|- 48\/A. 

Building on these lemmas, we can define the bond truncation procedure Truncate 
for our algorithm as follows. 

( 2 ) 

Truncate; Take each vector within S'/ ■, and truncate all bonds on 1,..., z — 1 
to P := 800nRi/8oon = Take S® to be 


:= I Truncp Truncp 


■ Trunc 


1-1 


s)||s) €Sg}uL;,,p 


As before, the following Claim establishes that the Truncate step preserves the via¬ 
bility with the advertised parameters. 

Claim 10 (Proof in Section 3.3). The set resulting from Truncate, is {i,pi 

q,P 2 ,d = l/5)-viable, where 

P 2 {'n) = n“ where a = 0(l . 


2.5 Error reduction 

Next step in the algorithm is to push the error level back down. To do this we need to 
construct operators called approximate ground state projectors (AGSPs). One convenient 
choice of AGSP would be of the form 


A 


:= exp 


x(B - eo)^' 
2e2 
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where x is a parameter defined later. We cannot however construct such an operator. 
Firstly it would require us to exactly know the ground state energy. Secondly it would 
require exponential time to construct and specify. 

We already have a ground state approximation I 71 ) which we got from the non¬ 
degenerate algorithm with energy at most eo+ 7 ^e/ 4 /. As such, if we define Cq := ( 7 i|i 7 | 7 i) 
then for sufficiently small error 7 we have 

|eo - e'ol < ^^e/4/ < e/2. 


In order to construct an approximation of A by a more efficient operator K, we define 

UB{rj), (4) 


. \T|r^ 

K := 4^ ^ exp 


'J2'KX 


j=0 




2x 


where Usit) is an approximation of exp{—iHt) with Schmidt rank B across every cut, the 
construction and analysis of which is given in Ref. [13]. 

The approximation between A and K is broken into three distinct approximations. 
The first is that the integral on (— 00 , 00 ) is truncated to an integral on [—T, T], with 
an exponentially small error in T. The second is that this integral is approximated by a 
Riemann sum, specifically the rectangle rule^ with a discrete step-size of r. The third is 
that exp(—iHt) is approximated by a unitary Usit) that has low Schmidt rank. All three 
errors are considered in Section 3.4. If we take parameters in Eq. 4 to have the scaling 

X :=0(logC"^) 

r:=0(e-Mog(l/C)v'log(n/C)) 

T~^ :=C>(n^C"Vlogn/C) 

R:=(C-i)‘^(iA).poly(n/C), 


then this approximate AGSP is capable of lowering the energy error down to A = As 
such we will henceforth refer to this as a /"-approximate AGSP. With the definition of K 
in hand, we can give a precise specification of Reduce: 

Reduce: Decompose a /-approximate AGSP with / = ce®//^ as K = Aj (g) 

Bj. Then return the set given by applying this approximate AGSP, combined 
with the recycled Schmidt vectors of the previous approximate ground states. 


c(4) ._ 


{ 


Aj |s) 


Vj, 



As we will prove in Section 3.4, the Reduce step can lower the error of a viable set 
from 5 = 1/5 down to A = /, increasing both the bond dimension and cardinality by a 
factor of 

(1//)^^^^/^) • {n/cr^^ . 

Therefore the above construction with / = ce^jf^ can send a {i,pi + q,P 2 ,S = l/5)-viable 
set to a {i,ppi + pq + q,pp 2 ,ce^/ set where 

p{n) = . 

This forms our next Glaim. 

Claim 11 (Proof in Section 3.4). The set resulting from Reduce, 44 ihPPi T pq + 
g,pp 2 ) A = ce^/-viable. 

^Higher-order approximations such as other Newton-Cotes formulae could be used, though this won’t 
change the overall scaling of the various parameters. 
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2.5.1 Final error reduction 


As for the final step, we now want to reduce the error further than was previously necessary. 
We do this by simply constructing a stronger version of the previously used approximate 
AGSP. 


FinalReduce: Construct a r7^/4/-approximate AGSP K. Return the set con¬ 
structed by applying this AGSP, combined with the previous approximate ground 
states. 



esS}u{|7i),---.l7»-i)). 


Claim 12 (Proof in Section 3.4). The set resulting from FinalReduce, is {n,pi + 
q -|- g,PoP 2 , A = rf / ^f)-viable, where 

Po{n) = n“, where a = 0[l + ydogrT^^Ylogn) . 


2.6 Ground state approximation 

Once we have a viable set on the entire system, we then need to extract the witness 
state. By construction, a viable set provides a domain on which known optimisation 
techniques can be applied to efficiently find approximate ground states. We summarise 
this with the last procedure. 

ApproxGroundState: Take cr to be a density matrix supported on Span{Sh,n), 
found as the solution to the convex program 

min Tr{Hah) (5) 

where i'yjlo'hhj) = 0 for 1 < j < h, 

Tr(<T,,) = 1, ah>0. 

Return \^h)i a leading eigenvector of a^- 

The objective functions and constraints can both be evaluated efficiently, firstly be¬ 
cause our witness states are MPS with small bond dimension, and the Hamiltonian is a 
matrix product operator (MPO) with small bond dimension, and also because the domain 
of the optimisation is polynomially large. Since the program takes the form of a convex op¬ 
timisation, these two properties are sufficient to allow the program to be solved efficiently. 
By optimality of the above program, the energy of Ch must be at most that of the witness, 
eo + rj'^efAf. If we take rj < 1/3 and using Lemma 15 we get that the leading eigenvector, 
| 7 /i), has an energy at most eo + rfe/^g, and ground space overlap ||G ly/i)!! > 1 — /^g- 


2.7 Proof of Theorem 1 

We can now turn to our proof of the main result Theorem 1. The proof is conditional on 
Glaims 6, 10, 11, and 12, whose proof is given below in Section 3. 

Proof. The Glaims 6, 10, 11, and 12 establish that in each iteration of Algorithm 1 a viable 
set with energy error at most ??^/4/ is produced, and Lemma 15 establishes that this yields 
an orthonormal set of vectors such that (yjj G {gj) > 1 — J^g for all j = 1 ,..., g. Define 
T := X]j=i l7j)(7jl! compute 

||G - Till = Tr(G - Tf = 2g - 2TV(GT) <2g- 2g{l - /2g) = g^ . 

This completes the proof of the main result conditional on the Glaims. □ 
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2.8 Runtime 

As all the parameters given in Table 1 are polynomials in n, and all operations required to 
execute Algorithm 1 (convex optimisation, inner product of MPS etc.) can be performed 
at polynomial overhead, the overall runtime is also polynomial in n. The leading order 
term in n gives a runtime of T = n“ where a = 0(l + y\ogrf~^^\ogT^, which reduces to 

J where rj = , 

where rj = . 

In Ref. [22], an improved analysis of the FinalReduce step (which we will discuss in 
Section 3.4.2) is presented. The scaling the polynomials po, pi, and p 2 are all improved to 
q = in the r]~^ = case. This means that the corresponding run-time reduces 

to T = 

One thing to note is that this leading order scaling is independent of the degeneracy. 
If we now focus on the degeneracy scaling, hxing n and allowing g to vary, we hnd that 
the dominant term now comes from the cardinality of the boundary contraction net. As 
such the leading-order scaling becomes 2^ where b = ; this is also the leading order 

scaling in the inverse-gap. 


3 Analysis 


3.1 Preliminaries 


To prove the claims in the previous section, we first state some simple results that are 
easy to prove. These results were also used in Refs. [1,2], and we state them here without 
proof. 

Lemma 13 (Overlap results [1,2]). 

• For any states for which |(u|rc)| > 1 — <5 and |(r;'|r(;)| > 1 — 5', then |(u|u')| > 

l-2(,5 + ,5')- 

• For any two states such that |(u|rc)| >1 — 5, |(u|0|u) — (r(;|0|r(;)| < 2\/^ for any 
operator O with ||0|| < 1. 

• If |(w|r)| < u) and |(u|tt)| > 1 — d then |(f |r)| < to + 


• For a pair of vectors with bounded overlap |(tii|n 2 )[ < uj, and a second pair [ui), |u 2 ) 
such that \{uj\vj)\ > 1 — <5 for j = 1,2, then we have |(ui|u 2 )| <u + pl05. 

Lemma 14 (Orthogonalising lemma). Suppose you have two vectors |tt) and |u), such that 
{u\H\u) , {v\H\v) < eo +Ae, with bounded overlap |(tt|u)| < (5. Then there exist orthogonal 
vectors \u '), \v') G Span{|n) , |u)}, such that {u'\H\u '), {v'\H\v') < cq + Ae(l + f3)/{l — fl). 


Proof. The proof is constructive, take the two vectors to be 

\v) — |u) {u\v) 


\u') = \u) , |u') = 


|u) — |li) {u\v) II 

As \v') oc (/ — |tf)(rt|) \v) this vector is clearly orthogonal to |tt'). As for the energy bound, 
dehne H' := H — cq, then 

{v\H'\v) — 2Re {v\H'\u) (rt|u) -|- (u|«) {u\H'\u) (rt|u) 


{v'\H'\v’) = 


|u) — \u) {u\v) 


< l±|Mhi At < 1±^ At 


1 — (w u) 


1-/3 
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After reintroducing the cq term in H we get the hnal energy bound. □ 

If we apply the same argument to the projector orthogonal to the ground space, the 
same (1 + /3)/(l — /I) growth can also be shown to occur in the overlap error. The Schmidt 
rank of the resulting orthogonalised vectors will in general be the sum of the two original 
vectors’ Schmidt ranks. If we can bound the overlap between a vector and any member 
of a vector space, then that vector can be similarly orthogonalised with respect to that 
whole space at the same cost. 

As expectation values are bilinear in pure states, the optimisations in the Trim and 
ApproxGroundState are both performed over mixed states such that their objective 
functions are linear. To aid in the analysis of these optimisations, we will next show that 
the leading eigenvector of a low energy mixed state must be a low energy pure state. 

Lemma 15 (Demixing). If a mixed state a has an energy at most cq + Ae/{2g-\- 1) where 
A < 1/3^, then the leading eigenvector has energy at most cq + Ae. 

Proof. First decompose the mixed state as cr = \vk){vk\j where Xk > 0. By the 

normalisation of cjfc we have ^ = 1. Next define the set of indices K as 

K = {k \ {vk\H\vk) < eg + Ae} . 

By comparing the energy bound of the mixed state as a whole as well as those on indices 
in A', we get 


eg + Ae/{2g + 1) > Tr(crif) 

= ^ Afc {vk\H\vk) + X] Afc {vk\H\vk) 
k&K k^K 

> eg + Ae Xk 
k(^K 

which in turn gives Ylk&K^k A ‘2‘gl{2,g + 1), analogous to a Markov inequality. As we 
required A < l/S^, each vector corresponding to a index in K has an energy at most 
eg + e/3g', and thus a ground space overlap at least 1 — 1/3^. By Lemma A.3 there can 
only be g such orthonormal vectors, leading us to conclude that \K\ < g. Using this we 
can bound the largest eigenvalue corresponding to an index in K 


maxAi. > meanAfc 

k£K k&K 


Y^keK Afc ^ ‘2, 

\K\ -25 + 1 


as well as that not in K 


max 

k^K 


Xk ^ 'y ^ Xk 

kf^K 


< 


1 

25 + 1 ■ 


As such we can conclude that the largest eigenvalue must belong within K, meaning that 
the leading eigenvector has energy at most eg + Ae. □ 


3.1.1 Non-degenerate ground state approximation 

Here we will briefly address NondegenerateGSA, which we have claimed is a small 
modification on the nondegenerate ground state approximation algorithms of Ref. [1,2]. 
The analysis of these algorithms largely holds in the degenerate case, appropriately re¬ 
placing the notions of ground state overlap with ground space overlap. The only key 
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difference is the use of a special case of Lemma 15. As Lemma 15 will be invoked during 
the size-trimming and ground state approximation steps (Sections 2 . 6 , 3.2), the energy 
levels in these steps grow by a factor of g overlooked in the non-degenerate analysis. By 
applying a slightly more powerful error reduction, dropping error levels beforehand by g 
to compensate, the same overall scaling can be proved, at the mere introduction of an 
extra constant factor. Due to this rather superficial change we will use these algorithms 
without proof of their validity. Furthermore we will take 

q{n) := re“ where a = 0(1 -|- e“^-\/logr/“^/logn) ( 6 ) 

to denote the bond dimension of this first approximate ground state I 71 ) given by this 
algorithm. We will show that the bond dimension of all subsequent ground state ap¬ 
proximations follows the same scaling, and thus use q{n) to generally denote the bond 
dimension of any combination thereof. 

3.2 Size trimming 

Next we will consider the size trimming sub-step, and prove the associated Claim 6 . In 
Lemma 19 we will show that there exists a low-energy state orthogonal to the existing 
ground states, which we hope to find. After this we will introduce and define local- 
distinguishability, and give some intuition about how our size-trimming works in the low- 
and high-distinguishability regimes. We will then consider how the error levels behave in 
these two regimes, showing that the size-trimming procedure yields a good viable set for 
any level of distinguishability. We will conclude this section in with a full proof of Claim 6 . 

3.2.1 Low-energy state 

The optimisation Eq. (1) in Trim is performed with respect to Hl and not the whole 
Hamiltonian. There is no reason to assume a priori that states which have low energy 
with respect to the global Hamiltonian H, have low energy with respect to half-chain 
Hamiltonians Hl or Hn. Using the locality of the Hamiltonian however we will see that a 
global ground state of H can be excited by an energy at most 1 with respect to Hl + Hr- 

Lemma 16. For any ground state |r) of H, 

(r| { h'l + h ' j ,) |r)<i 

Proof. Consider the state |u) := \l) ( 8 > |r), where \l) and \r) are in the 0-eigenspace of Hf 
and H'j^ respectively, which are non-empty by construction. The energy of \v) is 

{v\H\v) = €L + eR + {v\Hi\v) , 

so for |u) to not have an energy lower than that of a ground state, we require eq < 
Cl + er + {v\Hi\v). By bounding this further by the normalisation ||H^i|| < 1, this means 
the ground expectation of H'^ + H'^ is at most 1 . □ 

As global ground states have low energy with respect to Hl + Hr, there must exist a 
good approximant of the ground state in the low-energy subspace of Hl +Hr. As such we 
will restrict our analysis to this subspace, which we will refer to as the truncated Hilbert 
space. 
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Definition 17 (Truncated Hilbert Spaces). Let Pt be the projection onto ^], 

where „] is the subspace of spanned by eigenstates of H'j^ with eigenvalues 

at most t, that is, 

im(Pt) = span {|u)|H)j |u) = A |u) ; A < t} . 

Similarly let Qt be the projection onto (7^[i,i] <8) which corresponds to the span 

of the eigenvectors of + H'j^ with eigenvalues at most t, that is, 

im{Qt) = span{|u)|(H):j + iL^) |u) = A|u) ; A < t} . 

As the left and right Hamiltonians commute and are positive, we have Pt > Qt- 

The following choices for t and ^ ensure that the error of the boundary contraction 
step remains small. Define t := 99 (log 99 + 4 log / — log c — 6 log e) such that 99 • = 

gg6/j4 Purthermore without loss of generality we can take e < 10® < ^/TOf^Jc such that 

t > 1 + > {T\iHi + i7;j)|r) + 

for any ground state |r). Finally we define the spacing of the boundary contraction net ^ 

? := c'e/f{l + t) = n{{e/g)/log{g/e)). (7) 

In general whilst a ground state can be heavily frustrated, containing contributions 
from states with arbitrarily high local energy, these contributions will be small. Specifically 
we expect that a given ground state will have a large overlap with Pt and Qt — this result 
is known as the truncation lemma, and was proven in Refs. [2,12]. 

Lemma 18 (Truncation lemma [2,12]). For any ground state [F) of H, 

||(l-Pt)|F)||<||(l-gt)|F)||<99-2-*/®®. 

Given these results, we can now move on to proving the existence of a low-energy 
witness state which Eq. (1) will approximate. 

Lemma 19. There exists a state [tc) in Span{S'|^^}(8)'H^^j^ orthogonal to the previously 
constructed approximate ground states [ 71 ), ..., \'^h-i), of Schmidt rank B := Bs where 
5 = with an energy at most cq + c'elf. 

Proof. Let |u') be a witness of such that it is orthogonal to the previous ground state 
approximations and has energy at most eo + ce®//"^. By Lemma 3 there exists a ground 
state |F) such that |(u'|F)| > 1 — ce®//"^. Using Lemma 13 we get that 

{v'\Hi\v') > (F|H,|F) - 2V^e^/f . 

Using this together with the energy of \v'), and again assuming that ce^ff'^ is sufficiently 
small gives 

(u'l {Hf + H'j,) \v') < (F| [Hf + H'lt) |F) + 2V2ce^/f + ce^f 
<{T\{Hf + H'^)\T) + 4V~ceyf. 

Decomposing the above using the projector Pt and its orthogonal complement - both 
of which commute with Hf and - then this gives 

+ 77)j)|F) + + H'j,)\v') 

>{v’\Pt{Hf + H'j,)Pt\v') 

+ {v'\{l-Pt){Hf + H'j,){l-Pt)\v') 

> {v\{Hf + \\Pt |u')|f + 111(1 - Pt) 
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where we have defined |r;) := Pt \v') / \\Pt |wO||, such that |u) now lies in the desired domain 
SpanjS'l^^} ® by construction. Using ||Pt|u')|p + ||(1 — Pt)|u')|p = 1 as well as the 

inequality t > (r|(P£ + P^)|r) + Ay/ce^jf'^, we can reduce the above simply to 

+ H'r)\v) < (r|(Pi + p;,)|r) + 4v^eV/' . (8) 

The overlap between this truncated witness and the ground state is bounded 

|(^^|r)|>|(u'|Pi|r)| 

>Ku'|r)|-Ku'|(i-p,)|r>i 

> 1 - ce®//"^ - 99 • 

= l- 2 ce®//^ 

and so once again applying Lemma 13 gives 

{v\Hi\v) < (r|Pi|r) + . (9) 

Combining Equations ( 8 ) and (9) we can bound the total energy of the witness projected 
into the truncated Hilbert Space 

{v\H\v) < eo + ^yfce^/f . 

Next we need to considering trimming the bond dimension. By applying Lemma 9 and 
recalling the definition of from Eq. (3), we get that the truncated state 

|rc) := Trunc^ |u) / UTrunc^ |u)|| 

has {w\H\w) < eo + 96'y4ce// < eo + 0.75c'e// . 

Next we consider the overlap induced with the existing ground states. To bound this 
we consider the overlap between |t(;) and |u'). We had that |(u'|r)| > 1 — ce®//^ and 
|(u|r)| > 1 — 2ce®//^, Lemma 13 therefore gives 

|(u'|u)| > 1 — 6 ce®//"^. 

Next take I 7 ) := Trunc^ |r) / ||Trunc^ |r)||, which has |( 7 |r)| > 1— 8 -y/ce^//^ by Lemma 7. 
Lemma 13 once more gives that |(u| 7 )| > 1 — 20y/ce^jf'^ . As \w) is a truncation of |u), by 
the Eckart-Young theorem it must have a larger overlap than I 7 ) 

|(u|u;)| > |(u| 7 )| - 20 v^e^//^ . 

Applying Lemma 13 a third time we get finally that 

|(u'|w)| > 1 - . 

As I 71 ) ,..., | 7 /i-i) are all perpendicular to |u'), the overlap between |rc) and any element 
of Span{| 7 i) ,..., | 7 /i_i)} is upper bounded by the component of Ire) orthogonal to |u') 

y^l - |(U|w)P < \/l- ( 1 -52v^e3//^)^ < a/104v/^// . 

Taking e < 10^ without loss of generality this overlap is less than 1/7, so we can 
orthogonalise |u)) to the existing ground states at a multiplicative energy cost of 4/3, 
resulting in a final state with {w\H\w) < eo + e'e//, and all other desired properties by 
construction. 

□ 
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3.2.2 Distinguishability 


Viable sets capture local information about a witness state, specifically the Schmidt vec¬ 
tors. If two states are globally orthogonal, this can mean that their Schmidt vectors span 
orthogonal spaces, but they can also span identical regions. Indeed global orthogonality 
of witness states cannot be locally imposed for this reason. 

When considering viable sets, a natural notion of distinguishability is the overlap of 
the desired witness state, and the Schmidt vectors of the already constructed states. If 
we take \w) to be a witness state which we wish to capture, to be the set of Schmidt 
vectors on [l,i] of the h already constructed ground states, and Ilh,i the projector onto 
Span(L/j^j), then the distinguishability is dehned as 

Dh,i := 1 - {w\Uh^i\w) . 

From here we will omit the subscripts on L, 11, and D, since they will be clear from 
context. If we think of 11 as being the projector onto the vectors on [1, i] which we already 
have, then D can be considered a measure of how much of the state is unaccounted for. 

Clearly if D is low, then a viable set can be constructed entirely by recycling existing 
Schmidt vectors. The precise bounds for the quality of viable set constructed in this 
manner, as a function of D, are given in Lemma 20. 

In the case of a high D new vectors need to be found. In this case we know there is 
little overlap between the Schmidt vectors of our new witness and Span(L). Our proce¬ 
dure works by restricting to a low-energy subspace, and projecting away from the existing 
vectors, by minimising the expectation value of 11. In Lemma 21 we will show that this 
procedure must produce a low energy state. Given that our witness is highly distinguish¬ 
able, Lemma 22 will show that this new witness after size trimming must also be highly 
distinguishable. This will culminate in Lemma 23 where we show that this therefore im¬ 
plies that the new witness has small overlap with the previous ground states, producing a 
good viable set. 

We will conclude this subsection by showing in Lemma 24 that by performing both 
the low- and high-distinguishability procedures and combining the results, we can upper 
bound the error of our viable set in a distinguishability-independent manner. 

3.2.3 Low D 

By definition, in the low D case the left Schmidt vectors of |rc) have large overlap with L. 
As a result we can see that L will form a viable set, whose witness is the projection of Itc) 
onto L, and that the error of this viable set is low for low D. 

Lemma 20 (Low D Solution). There exists a state |u) orthogonal to I 71 ),..., | 7 h-i) with 
Tr[j_|_i „] |u)(u| supported on the span of L^ i, sueh that there exists a ground state |r) with 
|(r|u)| > 1 - (Ilow where 


^Low — 1 


1 - c'/f - y/D 


( 10 ) 


Proof. The idea here is to project Itc) into Span (Lh^i) and use the distinguisha¬ 

bility to bound the error associated with renormalisation. Specifically take 


v) 


n itc) 
^H\\ 
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where we recall that (t(;|n|r(;) = \ — D. As the existing ground state approximations are 
all stable under 11 by construction, the orthogonality of Ire) to these vectors carries over 
to je). The energy of |re) is at most eo + cT//, by Lemma 3 this means there is a ground 
state |r) with 

IHr)|>i-c7/. 

Decomposing this with the projector 11, and using both the triangle inequality and Cauchy’s 
ineqality, we get 

l-c7/<|(r|u;)| 

<|(r|n|u;)| + |(r|(/-n)|u;)| 

< I(r7)| + \/d. 

Rearranging this gives the desired ground state overlap stated above. □ 


3.2.4 High D 

Take a to both be the solution to Program 1 for the X closest to cont(tc) and the smallest Y 
which is greater than or equal to {w\Hi\w) — (yilLf^lyi), i.e. the tightest set of constraints 
satisfied by |rc)(u)|. By the nature of the two nets, described in Section 2.3.2, we have that 

||A — cont(r(;)||]^ < ^ and ^ < Y — {w\Hi,\w) + < Ac e/f . 

Here we want to argue that if D is not too low, then one of the a outputs of Eq. (1) 
has an eigenvector which can be extended to a good witness. To extend the eigenvectors 
of a we take the ‘boundary-uncontracted’ version a' = UwaUl,, eigendecomposing these 
matrices as: 

0 -= = '^^k\v'k){v'k\ 

k k 

where \v'^) = \vk)- 

In Trim we kept all eigenvectors whose value exceeded /g. We next want to show 
that there can only exist a finite number of such eigenvectors, and that they represent the 
vast majority of the spectrum of a. 

Lemma 21. There exist at most g eigenvectors of a' with energy at most eo + 30,OOOcT/5, 
and the sum of the corresponding eigenvalues is at least A = 1 — 1/(10,OOOgf + 5000). 

Proof. Dehning a' = U^jO'Uw, the boundary contraction property Lemma 5 gives an energy 
bound on a'. 


Tv{Ha') < TiiHLa) + {w\{Hi + Hr)\w) 

+ ||Tr[^j_i](cr) - cont('u;)||^ (^1 + UIHrU^ - {w\Hr\w)'^ 

< Tr(iLLCT) + (tcKiLi + iL/j)|r(;) + ||Tr[i^j_i]((T) - cont(r(;)||^ (^1 + uIHrU^ ) 

< {w\H\w) + {Y - {w\Hl\w) + (yiliLilyi)) + ^(1 + t) 

< eo + 6c'e// . 


By applying an argument similar to Lemma 15, if we increase the excitation not only by 
a factor of {2g + 1) but also by 5000, we get that there are between 1 and g eigenvectors 
which have energy at most eo + 30,000cT/5'. Moreover we also get that the sum of the 
corresponding eigenvalues is at least 1 — 1/(10,OOOgf + 5000) = A. □ 
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Next we rearrange the order of the eigenvectors. Let 1,... ,g refer to these g vectors^ 
described in Lemma 21. In other words: 


< Co + 30,000c'e/g' for A: < 5 
> Co + 30,000c'e/5( for A: > gr 

Next we drop any such vectors with eigenvalues less than the 10“^/^ threshold^. It should 
first be noted that by applying Lemma 15 to cr we get that the dominant eigenvalue that is 
at least 2/ {2g + l) has a corresponding eigenvector in this set, and so this set is non-empty. 
Next we note that as there are at most g — 1 such vectors dropped, so decrease in the sum 
of eigenvalues is at most 10“® • < 10“®. As such we can conclude that the sum of the 

corresponding eigenvalues is still bounded 

> A-10-9. 

k=l 

Without loss of generality, take |u^) to be the member of ,..., for which the 
expectation value of 11 is the lowest. We now want to show that the expectation of 11 with 
respect to |u^) not significantly higher than it is for Ire) (which is 1 — D by definition). As 
such this will imply that in the high D case |u() has small overlap with L, meaning that 
we have succeeded in finding a state with low overlap with the existing ground states. 


Lemma 22. 


W|n|ul) < 


T + 
A 


(1-^) 


Proof. The specific iteration of Eq. (5) we are referring to is that for which \w) satisfies 
all of the constraints of the optimisation. As such, by the optimality of the solution a, we 
know that (t(;|n|r(;) > Tr((Tn). Decomposing this and using the definition of D: 


(ui|n|u;) > Tr((Tn) 

1 - D > ^ Afc {vk\Il\vk) . 

k 


As n acts strictly to the left of the cut and strictly to the right they comute, meaning 
the expectation value of 11 for primed and unprimed eigenvectors is identical. 


l-D>Y,>^k{v'kW^) 

k 

9 

> ^Afc(4|n|4). 

k=l 

As we defined |u() to have minimal expectation this can be further loosened to 

^Padding out these indices with zeros as necessary if there are less than g such vectors. 

®Once again padding out with zero vectors such that there are still g such vectors for convenience. 
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Rearranging this together with the bound for Ylk=i we get 

(u'i|n|u'i> < [A - - D). 

After a binomial expansion, followed by a loosening of the inequality, we arrive at our final 
bound. □ 

In general |u^) will have some non-trivial overlap with the existing ground state ap¬ 
proximations. This low expectation value for 11 can however be used to bound this overlap 
when D is not too small. 


Lemma 23. For any vector I 7 ) G Span{| 7 i) ,..., | 7 /i_i)} the overlap between I 7 ) and |u^) 
is upper bounded 




T + 10”® 
A 


(l-D). 


Proof. The main property to use here is that IT I 7 ) = I 7 ) by definition. Using this and 
Cauchy-Schwartz gives 


|(7|ul)|' = K7|n|ul)|' < ||n = K|n|Ui) < 


A 


+ 10 


-9 


(l-D) 


□ 


Take \u) to be |u^) projected orthogonally to I 71 ) ,..., | 7 fe_i), as per Lemma 14. Then 
the ground space overlap of this vector can be bounded to be |(r|n)| > 1 — i^High for some 
ground state |r), where 


f^High = 


30,000c' 1 + /3 
~~9 W 


and 


/ 3 ^ < 


- + 10 
A 


-9 


(l-D) 


( 11 ) 


3.2.5 All D 

Next we can combine the two above errors. We have shown that there exist bounds on 
the overlap error 5 in the low- and high-D cases. By combining these, we can bound S 
independent of the unknown D. 

Lemma 24. The combined error 6 = min((lLow) f^High) over the total range of distinguisha- 
bilities is bounded 6 < 0.01 for all g >2 and D G [0,1]. 

Proof. Consider the case of D = 10“^ and g = 2. The two error levels in this case can be 
explicitly calculated, and are 

<5low( 5 = 2, D = lO"'^) Ri 0.00995050 < 0.01 
^High(5 = 2, D = lO-'^) Ri 0.00999947 < 0.01 

It can also be shown that for a hxed g that (5 low and Juigh are monotonically increasing 
and decreasing in D respectively, and that for a hxed D both decrease monotonically with 
g. As such we can conclude that 5 < 0.01 for all D G [0,1] and g >2. □ 

Claim 6 The set resulting from Trim, sf'\ is (i,pi, dsb + 5 = l/100)-viable, where 

7 'i(n), q{n) = where a = 0(1 + e~^\/ log 7 “^/log n). 
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Proof. As shown above, marginalising over distinguishability for degeneracies at least two, 
the overlap error is upper bounded by 1/100. The cardinality contribution comes in two 
parts: the results given by Program 1, and the previously obtained left-Schmidt vectors 
that are recycled. Let the cardinality of the Program 1 results be r where 

r= £1^ X 0{l/e) = 

E-values Energy net BC net Schmidt rank 

The cardinality contribution from recycled Schmidt vectors is bounded by the bond di¬ 
mension of each previous ground state approximation, as given in Section 3.1.1. This 
contribution is q{n) = n“ where a = 0(1 -|- y^logr/“^/ log n), and dominates over r. 
The overall cardinality is the summation of these two contributions, and as such is given 
by pi(n) = n“ where a has the scaling given above. 

The results of Program 1 are linear combinations of vectors in Span(S'|^^ U Lh,i). The 
contributions from both spaces to the bond dimension of the worst-case result is the prod¬ 
uct of their cardinalities and bond dimensions respectively. As such these contributions 
are dsb for and for L/j j respectively. □ 


3.3 Bond truncation 

The idea here is to use the existence of low-rank approximate ground states to argue that 
truncation induces little error due to its optimal nature as a low-rank approximation [16]. 
First we note an important result [1,2] about truncation, that large overlap with a low-rank 
vector is maintained under truncation. 


Lemma 25 ([1,2]). Given a vector [u) with Schmidt rank R across i\i -|- 1, for any \u) 

(ri-unP^/^iu) > |(u|u)| -6. 

(" 2 ) 

Let \u) denote a witness of 5) such that there exists a ground state [P) such that 


l(rk)| >1- 


100 


Fix |r) to denote the ground state with which |u) has maximal overlap. Next we invoke 
Lemma 7, taking r(re) ;= i?i/ 800 n such that there exists an MPS [ 7 ) with bond dimension 
bounded by r(n) such that 

|(F| 7 )| > 1- — . 

I\ I //I - gQQ 

Applying Lemma 13 we can combine these, showing [ 7 ) and |u) have large overlap 




/I 1 \ 18 

1-2 ^ + = 1 -^ 


V 800 100 


800 


Next take [re) to be the renormalised truncation of |u) across every bond 

Truncp ... Truncp“^ |m) 

||Truncp ... Truncp“^ |m)|| ’ 

where we have taken P{n) := 800nr(n). 
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Applying Lemma 25 across every bond we get that the overlap-error induced by trun¬ 
cation grows additively, giving a final overlap of 

19 

Kbk)! > l(7l«)l - n/(800n) > 1 - — . 

Once again applying Lemma 13 we get therefore that the ground space overlap of \w), 
which is lower bounded by the overlap with |r), can in turn be lower bounded 

|(r|u;)| > 1 - 2 ^— + —) = 1 -—. 

'' ' - VSOO 800y 20 

Once again we must turn our attention now to orthogonalising Ire) with respect to the 
existing ground states to satisfy the relevant orthogonality condition. The procedure for 
bounding this once again is to note that because |u) is orthogonal to Span{| 7 i),..., | 7 /i_i)}, 
any overlap |re) has with this space is upper bounded by the component of |u)) orthogonal 
to |rt). Applying Lemma 13 once again to their mutual overlap with |r), we get that 

,, , M /I 1 \ 6 

tekt) >1-2 -H- =1 - 

- VlOO 20y 50 

and thus that the orthogonal component is bounded 

^1 - |(u|u))|^ < v'l - 442/502 < 1/2. 

By Lemma 14 then we can orthogonalise jrc) with respect to the existing ground state 
approximations, at most tripling the overlap-error. Thus, after orthogonalisation, there 
must exist a ground state jT') 

|(r'|r(;)| > 1 --^> 1 --. 

I\ I /I - 20 “ 5 

Claim 10 The set resulting from Truncate, S') , is {i,Pi + q,P 2 ,S = l/5)-viable, where 

P 2 {n) = n“ where a = 0(1 -|- y^logr/“^/ log n). 

Proof. The error level 6 = 1/5 is proven above. As the truncation procedure is performed 
element-wise, the cardinality is changed only by the reintroduction of recycled Schmidt 
vectors, adding q. 

After truncation, the bond dimension is bounded by the trimming level 

P = 800n5i/8oon = . 

The orthogonalisation procedure then will additively increase this bond dimension by q, 
the bond dimension of the previous approximate ground states. As such the final bond 
dimension can be bounded by a polynomial p 2 of the same scaling. □ 

3.4 Error reduction 

Next we need to address the Reduce and FinalReduce procedures. We will do this 
by considering the construction of an AGSP A, and an approximation thereof K. Both 
operators are subject to a free parameter which controls the magnitude of the error 
reduction they can perform. The only difference between Reduce and FinalReduce 
will be the final level of ( utilised. As such we will analyse the procedure for an arbitrary 
C, stating the Reduce and FinalReduce as special cases of a more general analysis. 
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3.4.1 Exact AGSP 


First we consider the exact AGSP construction 


A := exp 


2e2 


where Cq := (7 i|i4|7i) is an inverse polynomial approximation to the ground state energy, 
and X is a ^-dependent parameter given in Appendix B. As well as specifying the precise 
parameters of the AGSP/AAGSP construction, in Appendix B we also show that for each 
ground state |r) and excited state IP"'"): 


l^|r)|| > 






It is in this sense that our AGSP approximates a ground state projector. 

Take \u) to be a witness of , such that there exists a ground state |r) with 


|(r|n)| >1-1/5. 


We can decompose this witness in the energy eigenbasis 

(i"-l 

1'“) = ’ 
j=0 

where IEq) is a ground state, {|£'j)}j>o are all excited states^ with energy ej > cq + e. 

(o\ 

Due to the role of \u) as a witness of we have po > 16/25, and due to normalisation 
YljPj — 1- Lastly take Ire) and to be the result of applying the AGSP to the witness and 
renormalising 

We can bound the renormalising denominator term by applying A to |n) in this de¬ 
composed form 


A|u)|| 


d"-l 


j=0 


> V^||A|Eo)|| - 


4 19 
“ 5 ' ^ 



'Y l-®i) 

i>i 


If we assume C < 1 then ||A|rt)|| > 1/5. Using this bound we can argue that because 
A shrinks the high energy terms, yet does not shrink the magnitude all that much, the 
energy of |tc) should be a drastic improvement over that of |tt). 


Lemma 26. 

{w\H\w) < eo + Ce/80,000 

'^As we only require that the vectors \Ej) be energy eigenstates up to rephasing, all vectors appearing 
in this decomposition from the same energy-eigenspace can be merged into a single vector, allowing for 
this decomposition to be considered non-degenerate without loss of generality. 
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Proof. We can bound the energy by simply expanding the expectation of H' := {H — cq), 
and seeing how it changes under the application of A using the above decomposition in 
the energy basis. 




{u\AH'A\u) 


<25 Y,Pj{EMH'A\Ek) . 

j,k>0 


Because A commutes with H, the double sum collapses and we have 


<25Y,Pj{EMH'A\E,) 

j>o 

< 25^pj(ej - eo)exp(-x(ej - ■ 

i>i 


Using Cq 


eo < e/2, and then using the Holder inequality, we have 

< 50^pj(ej - eo)exp(-x(ej - eo)Ve) 
i>i 


< 50 


max Eexpi—xE"^/e^) 
E>el2 


- 

^Pj 

-* 



The maximum is on the boundary because x is sufficiently large, so we find 

< 25ee"*/^ < Ce/3200 

Where above we have used the chosen value of x := 33 — 81og(C) from Appendix B. □ 


3.4.2 Approximate AGSP 


As well as the bounds on the power of A, Appendix B also gives that we can approximate 
A by a low bond dimension operator K with 


\\A-K\\<C,' 


Ce 

240,000n ■ 


( 12 ) 


For simplicity we can express K = A + C'O, where ||0|| < 1. As with K we can bound 
the renormalising term for |u), 


A|u)|| = ||(A + C'0)|u)|| 

>||A|u)||-C'||OWII 


Now if we take f < 1/10 then we get \\K |m)|| > 1/10. As such we can replace the state 
\w) that had an exact AGSP with an approximate version using the AAGSP, 


and similarly bound the energy. 


w') : = 


K\u) 

YWl 
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Lemma 27. 


(w'\H\w') < eo + Ce/400 

Proof. Once again we consider the expectation value of H' 

{w'\H'\w') ={w'\H’\w') 

_ {u\KH'K\u) 

\\K\u)f 

<100 {u\{A + C'0)H'{A+ C'0)\u) 

<100 {u\AH'A\u) + lOOC'^ {u\OH'0\u) 

+ lOOC' {u\OH'A\u) + lOOC' {u\AH'0\u) 

<4 {w\H\w) + 300nC^ 

<Ce/800 + Ce/800 = Ce/400 

Reintroducing cq gives the desired energy bound. □ 

Finally we must orthogonalise to negate any overlap with the existing vectors this has 
caused. 

Lemma 28. For any I 7 ) G Span{| 7 i) ,..., | 7 h - 1)}, |('«;'|7)l < 

Proof. The original witness |u) has a large ground state overlap, which in terms of the 
previous decomposition means 


y/m = {Eo\u) > 4/5. 


Using the fact that \\K |tt)|| < 1 and that A is diagonal in the energy basis, we can bound 
the overlap that \w') in turn has with \Eq). Using < 1/100, 


K^oh')l = 


\{Eo\K\u)\ 


..K\u)\\ 

>\{Eo\K\u)\ 


> 


VpJ {Eo\A\Ej 


-c' 


> {Eo\a\Eo) - c' 
4 19 1 _ 3 

“ 5^ ~ 100 “ 4 ■ 

Applying Lemma 13 we therefore get that 


l7K>|>i-2(i + i)=l, 


As has been argued before, the component of [re') overlapping with Span{| 7 i) ,..., | 7 /i_i)} 
can be bounded by the component perpendicular to \u). This component is bounded 



OVU ^ 400 -1 
10 “ 400 + 1 ■ 


□ 
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By Lemma 14 and the previous lemma, we can thus orthogonalise \w') to a vector \w") 
with a 400-fold increase in the excitation, resulting in an energy of 

< eo -|- ■ 

As such setting to the desired final energy error level and applying the local Schmidt 
operators of K, the energy error can be reduced to 

Claim 11,12 The set resulting from Reduce, is (i, ppi +pq + q, PP 2 , A = ce®//^)- 
viable, where 

p{n) = . 

The set resulting from FinalReduce, is (n, Pi + q + g, PqP2, A = 77^/4/)-viable, 
where 

Po{n) = n“ where a = 0(1 -|- y^log p~^/ log n). 

Proof. The construction described above gives a K with a bond dimension that scales (see 
Ref. [13]) as where o' = 0(1 -|- \/log C“^/ log n). 

For Reduce the cardinality growth is caused by the fact that K must be applied locally, 
i.e. each Schmidt operator must be individually applied. The number of these operators 
is bounded by the bond dimension of K. As such Reduce, for which ( = ce^ ff^ = 0(1), 
has a multiplicative growth in cardinality of p{n) = On top of this the recycled 

Schmidt vectors additively increase the cardinality once more by q. 

For FinalReduce however the whole operator can be applied without Schmidt de¬ 
composition, thus leaving cardinality unchanged. The recycled previous ground states 
grow the cardinality additively by g. 

The bond dimension increase is once again multiplicative, and once again bounded 
by the bond dimension of K as an MPO. For Reduce this means a growth by a factor 
of p{n). For FinalReduce however C, = i/^/4/ is not necessarily constant, as such the 
growth in bond dimension has the more general scaling of 

Pain) = n“ where a = 0(1 -|- y\ogr]~^^\ogn). 


□ 

Here we briefly comment on the improved analysis of Ref. [22]. Suppose we perform the 
following for our FinalReduce: instead applying of a single powerful AGSP, we interleave 
a mild AGSP with bond trimming multiple times. With sufficient parameters (the details 
are given in Sec.5.5 of [22]), this can be done such that the bond dimension remains at most 
'nP^^\ and the error level is halved. Performing this procedure the requisite number of 
times, the scaling of final bond dimension can be improved to po = for ri~^ = 

Garrying this analysis through, this also reduces the scaling of pi and q to both be 
and therefore a hnal run-time also of T = 

4 Conclusion 

We have shown the ground state approximation of Refs. [1,2] can be extended to degen¬ 
erate systems with an approximation of the ground space projector that is up to inverse 
polynomial error in any Schatten norm. 

AGSPs have proven powerful tools for relating the structural properties of gapped 
systems to various computational complexity bounds regarding them. Most AGSP con¬ 
structions are either functions of the Hamiltonians [2,12,14] or parts of the Hamiltonian 
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terms [1,17-20]. While we utilise such AGSPs in intermediate steps, our algorithm out¬ 
puts a state-based AGSP, and more thoroughly exploits the underlying one-dimensional 
structure, allowing for a more powerful overall construction. While previous AGSP con¬ 
structions have typically been full-rank, our AGSP takes the form of a minimal-rank 
projector, that moreover has a polynomially bounded bond dimension as a matrix prod¬ 
uct operator (MPO). It remains an open question for which systems there exist efficiently 
constructible AGSPs that are also low-rank AGSPs. Note that this is important for ap¬ 
plications because while AGSPs that are close to the ground space in operator norm are 
easy to construct even with small bond dimension in MPO form, they do not generally 
allow for low-error estimation of expectation values of other MPO observables. 

While running in polynomial-time, the specific runtime scaling of the above algorithm 
is highly suboptimal, leaving wide room for further optimisation. One potential avenue for 
optimisation is the level of redundancy within the viable sets that are generated. While 
there are rigorous upper bounds on the cardinality of these sets, no care has been taken 
to explore to possibility of high levels of linear dependence, opening the door to far better 
practical runtimes than our analysis would suggest. 

Note added in preparation: This work began as the undergraduate Honours thesis of 
the first author as a project to extend LVV to the degenerate case. During that time, the 
results of Huang in Ref. [2] (version 1) appeared, and were incorporated into our indepen¬ 
dent analysis. Near the completion of this manuscript, Huang also extended the results 
in Ref. [2] to the degenerate case and published this as Ref. [21] (version 3), an updated 
version of the original manuscript. While our techniques rely on the methods in version 1 
(Ref. [2]), our treatment of the degenerate case was made completely independently from 
version 3 (Ref. [21]). 

Our results and the new results in Ref. [21] are essentially the same, but with the 
following differences. The most notable is that the algorithm of Ref. [21] is single-pass, 
finding a full set of ground states with a single sweep of the system. In principle our multi¬ 
pass algorithm will generate viable sets with far higher levels of redundancy, causing worse 
scaling in the parameters of Table 1. In the final runtime however this slow-down only 
influences the various constants hidden under big-0 terms, giving an overall scaling that is 
identical to that of Ref. [21]. The second major difference regards the size-trimming, the 
most complicated step to generalise. Our algorithm works by applying an energy constraint 
and minimising the overlap between witnesses, while Ref. [21] constrains the overlap and 
minimises the energy. This allows one to use boundary contraction to circumvent the issue 
of local distinguishability, simplifying the analysis of this step. The final difference is that 
our algorithm does not require knowing a specific bound on the ground state degeneracy 
g, only that a bound of the form g = 0(1) exists, whereas it is not clear how to perform 
the algorithm in Ref. [21] without an explicit upper bound. 

Subsequent to the journal submission of this manuscript, but prior to its publication, 
Arad et.al. released a yet further improved algorithm in Ref. [23]. While our algorithm 
works by iterating along the spin chain, constructing matrix product states, the main new 
technique in Ref. [23] is to instead apply divide-and-conquer approach, generating states 
in a tree-like fashion. This method allows for much tighter runtime analysis (e.g. runtime 
reduction from poly(re) to 0(n^'^^) for frustration-free systems), as well as extensions to 
much larger classes of systems (e.g. a quasi-polynomial-time algorithm for certain gapless 
models). 

We thank Thomas Vidick for comments on the previous version of this manuscript. 
We acknowledge support from the Australian Research Gouncil via EQuS project number 
GE11001013, the US Army Research Office via grant numbers W911NF-14-1-0098 and 
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A Overlap lemma 

In this appendix we want to relax the notion of a basis, and show that key properties 
thereof are robust to a small amount of allowed error. 

Suppose we have some ^f-dimensional space, represented by a g-ia.uk projector G onto 
this space. Take \vi) for 1 < i < 5 to be an orthonormal set of vectors, each with 
{vi\G\vi) > 1 — 5/g. First we will prove that all members of S := Spanjjui)} have non¬ 
trivial overlap with G. We will then use this to show that the vectors G\vi) span the 
space. Finally we will show that any vector orthogonal to Spanjjui)} has low overlap with 
the space. 

Lemma A.l. Any vector ju) G S has (ujCju) > I — <5. 

Proof. As ju) lies in S, we can write it as a linear combination of the form 


9 



2=1 


Next we can bound the expectation value of the complementary projector I — G hy using 
the triangle and Cauchy-Schwartz inequalities. 


(u|(/ - Gjlu) = ^ c*cj {vi\{I - G)\vj) 


<'^\ci\\cj\\{vi\il - G)\vj)\ 



< S 
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where we have used the normalisation \ci\'^ = 1 and the resulting bound |cj| < y/d- 
This in turn gives {v\G\v) > 1 — <5. □ 

Lemma A.2 (Basis under projection). The projections G\vi) form a basis o/im(G). 

Proof. As the number of vectors matches the dimension of the target space, proving linear 
independence is sufficient to prove the formation of a basis. If the projected vectors were 
linearly dependent then this would imply there exist a non-zero |u) E S' such that G\v) =0. 
Lemma A.l however gives (u|G|u) > 1 — 5 > 0 for all such vectors, thus proving linear 
independence by contradiction. □ 

Next we want to generalise the concept of a full basis, specifically the fact that any 
vector perpendicular to a full basis is necessarily perpendicular to the entire space, meaning 
that there is a natural maximum to the size of any given basis. 

Lemma A.3 (Fullness). Any vector |u') which is orthogonal to S has {v'\G\v') < 5. 

Proof. As the projection of this space is spanning, we can pick |u) E S such that it projects 
into the same vector as |u'), i.e. G \ v) oc G \ v'). As such we can decompose both states as 

\v) = y/X \g) + Vl - A \h) 

1 ^;') = ^\g) + y/fTTxi\h') 

where \g) E im(G) and \h) , \h') E ker(G). The orthogonality of |u) and |u') gives 

0 = = VxVy + Vl — AVI — A' {h\h''j . 

Using \{h\h')\ < 1 in turn gives A + A' < 1. Lemma A.l gives that A > 1 — <5, leading us 
to conclude that A' < 6, and thus that {v'\G\v') < 5 generically for any vector orthogonal 
to 5. □ 


B Approximate AGSP construction 


In this section we explicitly lay out the various errors associated with approximating the 
AGSP defined in Sections 2.5 and 3.4. The scaling of these errors was stated in less detail 
in Ref. [ 2 ]. First we consider the AGSP itself 


A 


exp 


2e2 




r 

2x 


i{H 



dt 


where c'q is the ground energy approximation given by I 71 ). If we take < 1/2 and 
assume that Q > rf /df then 


rP'e 

|6'o-eo|<^<eC< 


e 

2 ■ 


The magnitude of this AGSP on the excited subspace, 
be bounded 



< exp 


fie/m 

2e2 


known as the shrinking factor, can 

= . 
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If we define x := 33 — 81og(^, then this shrinking factor is at most (^/2. Further assuming 
that ( < 1/25, the magnitude of A on the ground space can be lower bounded 


A\r)\\ =exp 
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> — . 
- 20 


3;(eo - eo)^ 

2 e 2 


Next we move on to the approximate AGSP, defined by 


K := 



\TM 

j=0 




2 x 


UniTj) ■ 


where T, r and D are parameters and is some unitary operator all of which we 
shall define below. This definition relies on taking the Fourier decomposition of A, and 
approximating this by truncating and discretising the integral and approximating the 
propagators by a known algorithm [13] . 


B.l Time Truncation Error 


The truncation error is due to the finite range of the integral, controlled by the parameter 
T. If we let 5t denote the truncation error (in the operator norm) then it can be bounded 


6t = 2 


V2 


TTX Jt 


exp 




2x 


— iHt 


dt 


< 


2 e 


y/2TTX 
e 


erfc 


J 


exp 




ff' 

2x 


dt 


< exp 


2 x 


where the last inequality holds for sufficiently large T, e.g. T > 1. For the time-scale 

T := ^ = 0(e-i log(l/C) v'log(n/C)) 

we can in turn upper bound 6t by 6t < CV^ = ll(C/^)j where (' is defined in Eq. (12). 


B.2 Discretisation Error 


Next we consider the discretisation error associated with approximating this integral with 
a Riemann sum. The term being approximately integrated is a Gaussian centred on the 
origin and so its absolute value is even as a function of t, allowing us to consider only 
the t > 0 error and doubling it. On this range the absolute value is not only even, but 
monotonically decreasing. As such the concepts of left and right Riemann sums correspond 
to upper and lower Riemann sums. If we let the discretisation error he 6 d then 





\TM 

^ rexp 
j=0 


2x 


iHrj 



r 

2x 


-iHt 


dt 
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This error can be upper bounded in turn by the difference between upper and lower 
Riemann sums which is equal to the left and right Riemann sums owing to the monotonicity 
of the integrand. 


< 2 


er 


^/2'kx 


\TM 

3=0 


2x 


— iHrj 


\T/T^+l 

Y """p 

i=i 


2x 


— iHrj 


The difference between these two summations is telescopic. As such it equals the difference 
between the highest and lowest terms. 


(5d < 2 


er 


V2 


TTX 


1 — exp 


2x 


-iHT 


Using the trivial bound on the real exponential term, we get 


< 2 


< 4 


er 


\l2'KX 

er 


1 — exp [ — iHT] 




TTX 


Again this error can also be bounded by 6d < CV^ if we take the time-step 

12 e 


T ^ : = 


12 /log(3/C') 


= O 


n 


C' V vr yC\/iogi/c 


B.3 Bond Truncation Error 

To upper bound the total truncation error for the entire approximation by CV^j if is 
sufficient to bound the truncation error for each propagator exp(—by since the 
total sum is normalized. Using the construction of Ref. [13] this error can be achieved if 
the bond dimension of each propagator is 

D = exp [O (rj) + O (log(nVCO)] 

< exp [O (T) -b O (logn/C)] 

= 2‘^(^ipoly(n/C) 

^ 2o(6-ilog(i/C)Vlog(VC)) poly(n/C). 

As bond dimension is additive, and there are 2T/(t -|- 1) terms, the total bond dimension 
for the AGSP B is 


B < 2Tt-^D 

= 0{e-^ log(l/C)Vlog(n/C))0 

^ 2o(6-Uog(i/c)Viog(VC))po|y(^/^) 

— Vlog(l/C)/ log(n)) 


^0{e Uog(l/C)Vlog(VC))po|y(^/^) 


where we have assumed C ^ 
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C Frustration bound 


Consider a ID system with nearest neighbour interactions, were each local Hamiltonian 
term Hi is bounded 0 < < 1. For a cut i\i + 1 we split the Hamiltonian into a 

left-Hamiltonian Hl = middle-Hamiltonian Hm = Hi, and right-Hamiltonian 

Hr = E •=.+! Hj] we refer to their expectation values as the left/middle/right-energy. 

Lemma C.l. For a Hamiltonian normalised as above, the left-energies relative to any cut 
for any two states with an energy less than AE over the ground, differ by no more than 
1 + AE. 

Proof Suppose we have two states |ui) and |u 2 ) with local-energies given in Table 2. 
where E = Eq AE, and Eq is the ground energy, and we take AL > 0 without loss of 


State 

(H) 

{Hl) 

{Hm) 

(Hr) 

bi) 

E 

L 

E-L-R-AR 

R + AR 

\V2) 

E 

L + AL 

E-L-R-AL 

R 


Table 2; Local energies of the two ground states across a bipartition, 
generality. Using the positivity of Hm on the expectations for \v 2 ) gives 

L + R< E - AL. 

As |ui) has a left-energy L, it must possess at least one left Schmidt vector |/) with a 
left-energy at most L, similarly \v 2 ) has a right Schmidt vector |r) with right-energy at 
most R. As such we can construct a state |u) = |/) (8) |r) with energy upper bounded 

{v\H\v) = {v\Hl\v) -P {v\Hm\v) + {v\Hr\v) 

= {1\Hl\1) + {v\Hm\v) + {r\Hji\r) 

<L+1+R 
< E + 1 - AL . 

As {v\H\ v) P Eq, we get that in general the difference between the left-energies is bounded 

AL <1 + AE. 

A simple example of a perturbed Ising model shows that this bound is tight. □ 
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